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Abstract 

We have studied the spin polarized hot neutron matter in the presence of strong magnetic 
Q ■ field. In this work, using the lowest order constrained variational method at finite temperature 

■ and employing ^IVig nuclear potential, some thermodynamic properties of spin polarized neutron 

matter such as spin polarization parameter, free energy, equation of state and effective mass have 

j. ■ been calculated. It has been shown that the strong magnetic field breaks the symmetry of the 

OO- 

■ free energy, leading to a magnetized equilibrium state. We have found that the equation of state 

o; 

. becomes stiffer by increasing both magnetic field and temperature. The magnetic field dependence 

I of effective mass for the spin-up and spin-down neutrons has been investigated. 
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I. INTRODUCTION 



Based on the supernova models, after the gravitational collapse of a degenerate stellar 
core and the explosive ejection of outer layers material, a protoneutron star would be born 



l|. Through the formation of the protoneutron star, the system arrives at a temperature 



Q 

of about 20 — 50 MeV [2]. This protoneutron star is hot, opaque to neutrinos, and larger 
than an ordinary neutron star [1]. After the formation of the protoneutron star, neutrino 
emission is the dominant process in cooling of the neutron star (mainly by URCA process 
and neutrino Bremsstrahlung) I2I. 

Woltjer has predicted a magnetic field strength of order 10^^ G for neutron stars as a 
result of the magnetic flux conservation from the progenitor star jsl- This in agreement 
with the experimental indication that the surface magnetic field strength of magnetars can 

^magnetar 



be of the order Bmagnetar ~ 10^^ - 10^^ G The magnetic field can be distorted or 



amplified by some mixture of convection, differential rotation, and magnetic instabilities 
[gI, I?]. The relative importance of these ingredients depends on the initial field strength 
and the rotation rate of the star. For both convection and differential rotation, the field 
and its supporting currents are not likely to be confined to the solid crust of the star but 
instead distributed in most of the stellar interior, which is mostly a fluid mixture of neutrons, 
protons, electrons, and other more exotic particles [sl- Thompson et al. argued that newborn 
neutron stars probably combine vigorous convection and differential rotation, making it likely 
that a dynamo process might operate in them j^. They expected fields up to 10^^ — 10^^ G in 
neutron stars with few-millisecond initial periods. In the core of high density inhomogeneous 
gravitationally bound neutron stars, the magnetic field strength can be as large as 10^" G 



10l |. In addition, considering the formation of a quark core in t 
neutron star, the maximum field reaches up to about 10^*^ G 10, 



le high density interior of a 



11| . According to the scalar 



virial theorem based on Newtonian gravity, the magnetic field strength is allowed to be up 
to 10^® G in the interior of a magnetar [12]. Moreover, general relativity predicts the allowed 
maximum value of the neutron star magnetic field to be 10^^ — 10^° G [13]. By comparing 
the cooling curves of neutron stars with the observational data. Yuan et al. obtained the 
magnetic field strength of order 10^^ G for many not so old neutron stars 14|. 

The finite temperature and strong magnetic field in the interior of a protoneutron star 
can influence different astrophysical quantities. Therefore, to have a better understanding of 
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different astrophysical phenomena such as supernova explosion, thermal evolution and cool- 
ing of protoneutron stars, gravitational wave emission spectrum from neutron star mergers, 
and to get the more precise astrophysical quantities such as properties of very young hot 
neutron stars and composition of neutron stars, one should consider the neutron star matter 
at finite temperatures (excited neutron star matter) and strong magnetic fields. 

Since /9-equilibrium leads to an increase in the number of neutrons in neutron star matter, 
it is possible to approximate the neutron star matter by the pure neutron matter. Many 



works have dealt with the study of dense neutron matter at finite temperature 15-|20|. 
Alonso et al. have used a field theoretical model for the analysis of relativistic neutron 



matter 



15| . By solving the model in the renormalized Hartree approximation, they have 



bund that an increase in the temperature moderately affects the 
l?! . In the Hartree- Fock approximation using Skyrme type interac- 



investigated the effect of central temperature on the maxima of mass for stable configu- 
rations, the radii of the configurations, and the gravitational red shift at the surface of a 
neutron star. Panda et al. used a mean-field description of nonoverlaping nucleon bags 
bound by the self-consistent exchange of a, u, and p mesons to investigate the properties of 
neutron matter at finite temperature [l^. They showed that by increasing the temperature, 
the equation of state becomes stiffer. Within the framework of the Brueckner-Hartree-Fock 
formalism and using the AVis nucleon-nucleon interaction for the spin polarized neutron 
matter, Bombaci et al. : 
single-particle potentials 
tions for spin polarized neutron matter, Rios et al. showed that the critical density at which 
the ferromagnetism takes place, decreases by temperature isj. Within the self-consistent 
Green's-function approach applying the CD Bonn and the AVis potential for neutron mat- 
ter, Rios et al. found that the effect of dynamical correlations on the macroscopic properties 
is rather insensitive to the thermal effects . The variational theory for fermions at finite 
temperature and high density has been applied by Mukherjee to neutron matter [20]. It has 
been found that the temperature dependence of the correlation operator is weak, but it is 
not negligible. Besides, it has been shown that the first order phase transition due to neutral 
pion condensation has a critical temperature of about 22MeV for neutron matter. The effect 
of strong magnetic field on the properties of dense neutron matter has also been considered. 
Perez-Garcia showed that for the neutron matter in the presence of strong magnetic fields, 
in the Skyrme model, there is a ferromagnetic phase transition at p ~ 4po, whereas it is 
forbidden in the DIP model 2l|. The results indicate that the effects of temperature on 
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the neutron magnetization remain moderate at temperatures up to about T = 40 MeV. In 
the context of the Landau theory of normal Fermi hquids, using Skyrme and Gogny effective 
interactions, some thermodynamical quantities such as isothermal compressibility and spin 
susceptibility of pure neutron matter have also been studied 

In our previous works, we have investigated the spin polarized neutron matter 23(|, sym- 



metric nuclear matter 



24( 1 . asymmetric nuclear matter as well as neutron star matter 25 1 



and magnetized neutron matter 



26| at zero temperature using lowest order constrained vari- 



ational (LOCV) method with the realistic strong interactions. We have also investigated the 

n 

thermodynamic properties of spin polarized neutron matter |27|, symmetric nuclear matter 
and asymmetric nuclear matter 29| at finite temperature with no magnetic field. In 



the present work, we calculate different thermodynamic properties of spin polarized neutron 
matter at finite temperature in the presence of strong magnetic field using LOCV technique 
employing AVis potential. 



II. LOCV FORMALISM FOR THE SPIN POLARIZED HOT NEUTRON MAT- 
TER IN THE PRESENCE OF MAGNETIC FIELD 

We consider a homogeneous system of interacting particles with A^*-"*"-* spin-up and 
iv(-) spin-down neutrons under the influence of a uniform magnetic field, i.e. B = Bk. 
The number densities of spin-up and spin-down neutrons are presented by p^~^^ and p*^^-* 
respectively. We introduce the spin polarization parameter, 6, by 

ni+) _ pi-) 

where — 1 < 5 < 1, and p is the total number density of system. The magnetization density 
of neutron matter is defined as 

m = (2) 

where pn is the neutron magnetic moment. The total magnetization of a given volume is 
also as follows 

M = ! mdV. (3) 
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In order to calculate the energy of this system, we use LOCV method as follows: we 
consider a trial many-body wave function of the form 

^ = ^0, (4) 

where is the uncorrelated ground-state wave function of independent neutrons, and T is 
a proper A^-body correlation function. Using Jastrow approximation [30] , T can be replaced 
by 

:F = S^j[i3). (5) 

where S" is a symmetrizing operator. We consider a cluster expansion of the energy functional 
up to the two-body term, 



ME 



The one-body term, is given by 

*=+ - 

where the first term of Eq. ([7]), shows the external magnetic field. It should be noted 
that we have used B to present the magnetic field strength; while the total magnetic field 
is the sum of the external magnetic field and the induced magnetization, B = H + AttM. 
Because of the tiny value of the neutron magnetic moment, we assume that the induced 
magnetization has a small contribution to the total magnetic field. Consequently, 

B ^ H. (8) 

Using Eqs. ([2]), ([3]) and (IE]), the one-body term can be written as 

^1 = -fXnB5 + J2 ^i- (9) 

i=+- 

The second term in Eq. ([7]), is as follows 

e. = J2^^^if'^T,B,p^^^), (10) 

k 

where nj(/c, T, B, p«) is the Fermi-Dirac distribution function in the presence of magnetic 
field, 

ni{k, T, B, pW) = ^^[,^(fc_r,s,pW)-7i,(T,i?,p(0)] (11) 

5 



In Eq. (fTTj) . li and JI^ are the single-particle energy of a neutron and the neutron chemical 
potential respectively. The single-particle energy, e,, of a neutron with momentum k and 
spin projection i in the presence of magnetic field is approximately written in terms of the 
effective mass as follows 13] 



2m;(T,p) 



+, 



(12) 



+ fXnB + U4T, p 



2ml{T,p) 

In fact, we use a quadratic approximation for the single particle potential incorporated in 
the single particle energy as a momentum independent effective mass. [/,(T,p«) is the 
momentum independent single particle potential. The effective mass, m*, is determined 



variationally 
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-|35|. The chemical potential, /ij, is also obtained by applying the constraint 



^n,(fc,T,i?,p«)=iV«. 



(13) 



The two-body energy, E2, is 



^2 = ^J2(^J W{12)\ 



(14) 



where 



K12) = -|^[/(12), [VL, /(12)]] + /(12)F(12)/(12). 



/(12) and V{12) are the two-body correlation function and nuclear potential respectively, 
n the LOCV formalism, the two-body correlation function, /(12), is considered as follows 



36|, 



/(12) = X]/^'Hn2)P, 



(fc) 
12 5 



k=l 



where 



3(fc=l-3) 

12 



'11 N /I 1 lo N .1 1 Ion 

- -^1.^2), (- + -ai.as + -S,2), (- + -^1.^2 - -5l2) 



(15) 



(16) 



In Eq. (fT6l) . 6*12 and cri.cr2 are the tensor and Pauli operators respectively. Using the 
above two-body correlation function and the AVis two-body potential [37| , after doing some 
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algebra, we find the following equation for the two-body energy: 

(2J+1) 



ttV (2m) ^^^^2(2^+1) 



[1 - (-1) 



L+S+li 



2'^zi2^z2 I SSz 



X 



/ 



X dr 



ZjTYI 9 

- + V; - 3V;, + 2(Vt - 3V;t) - 2V;ja« (r, p, T) 



fe=2,3 



2m 



+ + \4 + K + + (-6A; + 14)(\4, + Vt) -{k- l){Vlsr + 

+2[Vt + V,T + {-Qk + U)VtT - V;,]}a(f)'(r, p, T) 



where a = {J, L, S, Sz} and the coefficients a^'* , 6^, Ca^ , and are defined as 

a«'(r,p,r) = r^7i,5.(r,p,r), 



(17) 



(18) 



a(f)^(r, p, T) = r^[77^_i,5.(r, p, T) + ^/j+i,5.(r, p, T)], 



(19) 
(20) 



6^(r, p, T) = r2[/323/j-i,s.(r, p, T) - /323/j+i,5. (r, p, T)], 



cW^(r,p,r)=rVi7^,5.(r,p,r), 



cL^^'('^,P,2") = r^[?72/j-i,s.(r,p,T) + i/2/j+i,s,(r,p,T)], 



c(f)^(r,p,r) = r^[7^3/j-iA(r,p,r) + i/3/j+iA(r, p, T)], 



4^''(r,p,T) = r2[e2/j-iA(r,p,T) + A2/j+iA(r, p, T)], 



(21) 
(22) 
(23) 
(24) 
(25) 
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df\r,p,T) = r'[^sIj-i,sAr,P,T) + A3/j+i,5.(r, p, T)], 



(26) 



with 



^ 2J + 1' ^ 2J + 1' 



/3: 



23 



2J(J+ 1) 

2J + 1 



J2(J+1) J3 + 2J2 + 3J + 2 



2J + 1 



2J+ 1 



^2 



J(J2 + 2J+ 1) 
2J+ 1 



^3 



J( + J + 2) 
2J+ 1 



J3_^2J2 + 2J+1 J(j2 + J + 4) 

S2 — r, T , 1 ' S3 — 



2J+ 1 



2J+ 1 



_ J( J2 + J + 1) _ J3 + 2 J2 + 5 J + 4 

A2 — „ T , 1 ' -^3 — — 



2J+ 1 



2J+ 1 



(27) 



(28) 



(29) 



(30) 



(31) 



and 



{r, P, T) 



1 



27r6p2 



dfci rffca ni{k^,T,B,p^'^)nj{k^,T,B,p^^'^)J%\k2 - k,\r)- (32) 



In the above equation, Jj{x) is the Bessel function. 

Now, we minimize the two-body energy with respect to the variations in the function /< 
subject to the normahzation constraint. 



ii) 

a 



(33) 



The minimization subject to the above normahzation constraint leads to the normahzation 



of the two body wave function to unity 
function hs^{r) is defined as follows. 



36|. For the spin polarized hot neutron matter, the 



P 



-1/2 



where 



; 5, = ±1, 

; ^. = 0, 

dkni{km,T,B,p^^)Jo{kr)k''. 



(34) 



(35) 
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From the minimization of the two-body cluster energy, we get a set of coupled and uncou- 
pled differential equations js^. By solving these equations, we can obtain the correlation 
functions to compute the two-body energy term, E2. As the final step, we calculate the 
free energy per particle, F, to get different thermodynamic properties of spin polarized hot 
neutron matter, 

F(p, T, B) = E{p, T, B) - TS{p, T, B), (36) 
where S is the entropy per particle, 
S{p,T,B) = E 

«=+,- k 

+ n,(fc,T,5,p«)ln[n,(A;,T,5,pW)]}. (37) 



It should be noted that in our calculations, we introduce the effective masses, m*, as vari- 



ational parameters 3lN35|. We minimize the free energy with respect to the variations in 
the effective masses, and then we obtain the chemical potentials and the effective masses 
of spin-up and spin-down neutrons at the equilibrium state. In our approach, the effective 
mass depends on both density and temperature but it is independent of the momentum. 
The effective mass of a quasiparticle near the Fermi surface for the spin polarized neutron 
matter at low temperatures is also the static physical quantity of interest in the context of 
Landau Fermi liquid theory 



III. RESULTS AND DISCUSSION 



In Fig. [H we present the free energy per particle of spin polarized neutron matter 
versus the spin polarization parameter S. Fig. [1^ shows that in the presence of magnetic 
field, the free energy is not a symmetric function of spin polarization parameter and the 
equilibrium configuration would experience a net magnetization. Clearly, the effects of 
magnetic fields below B ~ 10^® G are almost insignificant, but by increasing the magnetic 
field, the equilibrium value of the spin polarization parameter (i.e. the polarization that 
minimizes the free energy) and the free energy decrease, leading to a more stable system. 
Fig. [T]d indicates that the effect of temperature on the fully spin polarized neutron matter 
is less than that of the unpolarized one. 
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Fig. [2] presents the equilibrium value of the spin polarization parameter versus density p. 
Fig. [2^ shows that at low densities (p < 0.2 fm"^), the magnitude of the spin polarization 
parameter decreases by increasing the temperature. However, at higher densities, the related 
values of the spin polarization parameter at different finite temperatures are almost identical 
to the one for zero temperature. This is due to smaller values of T/e*j- at high densities (Fig. 
Wp), in which e*j is defined as follows, 

-f-Y. (38) 



with 



4+ = - /^"^' (39) 



and 



.^_ = ^ + p„5. (40) 

In the above equations, e*^^ and k^p = (Qir'^p^^^)^ are the Fermi energy and Fermi momentum 
of neutrons with spin projection i in the presence of the magnetic field. It is evident from Eq. 
( l38l) that e*j gives an average of the Fermi energy of magnetized neutron matter. Therefore, 
the ratio T/e*j is a criterion for the fraction of particles which are thermally excited 40|, and 
how much the system is disordered. In Fig. [3l we show the spin polarization parameter at the 
equilibrium state as a function of the magnetic field B. At each temperature, the magnitude 
of spin polarization parameter grows by increasing the magnetic field (Fig. [3^). We have 
found that at strong magnetic fields, the effect of finite temperature is more significant 
because the ratio T/6*jr rises with the increase in the magnetic field (Fig. |3b). 

The free energy per particle at the equilibrium value of the spin polarization parameter 
is presented in Fig. |H It can be seen that at finite temperature, the free energy is an 
increasing function of density (Fig. H^). At low densities, the rate of increase in the free 
energy varies by increasing the density, but at high densities, this rate of increase is nearly 
constant. Moreover, the effect of temperature on the free energy is more pronounced at low 
densities. Fig. Hb shows that the free energy decreases by growing the temperature nearly 
at the same rate for different magnetic fields. The free energy decreases by increasing the 
magnetic field (Fig. Ht). We can see that by increasing the magnetic field up to a value of 
about B ^ 10^^ G, the free energy per particle slowly decreases, and then it rapidly decreases 
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for the magnetic fields greater than this value. This indicates that, above B ~ 10^^ G, the 
effect of magnetic field on the free energy of the spin polarized neutron matter becomes more 
important. 

From the free energy per particle of magnetized neutron matter, F, we can obtain the 
corresponding pressure of neutron matter using the following relation. 

This equation of state (EoS) is plotted in Fig. |5l For each value of the density, pressure 
increases by growing the magnetic field (Fig. [5^). This stiffening of the equation of state is 
due to the inclusion of neutron anomalous magnetic moments. From the astrophysical point 
of view, it should be noted that this stiffening of the EoS leads to the larger value for the 
maximum mass of neutron star 



4l|-|43l]. At each density, the pressure at finite temperature 
is larger than that of zero temperature (Fig. \5jp). It means that the equation of state of 
neutron matter becomes stiffer by increasing the temperature. Fig. ^ also shows that 
by increasing the temperature, the pressure increases nearly at the same rate for different 
magnetic fields. 

Fig. E] shows the effective mass corresponding with the equilibrium of the system for the 
spin-up and spin-down neutrons as a function of the magnetic field B. At low magnetic 
fields, the effective masses of spin-up and spin-down neutrons are nearly identical because 
the effective masses of spin-up and spin-down neutrons have the same values at 5 ^ 0. Fig. 
|6] indicates that the effective mass of spin-up (spin-down) neutrons decreases (increases) by 
increasing the magnetic field in agreement with the results obtained in Ref. [221 . From the 
comparison of Fig. |3K and El we can see that the shift in mass is due to the polarization 
of neutron matter. For the maximum value of the magnetic field considered in this work, 
i.e. 5 X 10^* G, and at T = 10 MeV and p = 0.3 fm^'^, the equilibrium value of the spin 
polarization parameter is about Sb^^^ = —0.23 and that corresponds to a mass shift of an 
amount of A(m*/m) ^ 0.02 with respect to the unpolarized case. 



IV. SUMMARY AND CONCLUDING REMARKS 

We have investigated the effect of strong magnetic fields on the thermodynamic properties 
of spin polarized hot neutron matter applying LOCV method and using AVis potential. We 
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have found that in the presence of a strong magnetic field, the free energy is not a symmetric 
function of the spin polarization parameter and the system is macroscopically magnetized. 
By increasing both density and temperature, the magnitude of the equilibrium value of the 
spin polarization parameter decreases. At low magnetic fields, the free energy decreases 
very slowly by increasing the magnetic field, but at stronger magnetic fields, the free energy 
decreases rapidly with the increase in the magnetic field. It has been found that the equation 
of state becomes stiffer by increasing the magnetic field. This stiffening of the EoS leads to 
the larger value for the maximum mass of neutron star. 
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FIG. 1: Free energy per particle as a function of the spin polarization parameter S: (a) for the 
cases B = G (solid curve), B = 10^^ G (dashed curve) and B = 5 x 10^^ G (dashdot curve) at 
the fixed values of the temperature, T = 10 MeV, and the density, p = 0.2 fm~^, (b) for the cases 
T = MeV (sohd curve), T = 10 MeV (dashed curve) and T = 20 MeV (dashdot curve) at the 
fixed values of the magnetic field, B = h x 10^^ G, and the density, p = 0.2 fmT^. 
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P (fm ') 



FIG. 2: (a) Spin polarization parameter at the equilibrium state as a function of the density p 
for the cases T = MeV (solid curve), T = 10 MeV (dashed curve) and T = 20 AleV (dashdot 
curve), and a fixed value of the magnetic field, B = 5 x 10^® G. (b) Same as in the top panel but 
for the ratio T/e*.. 
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FIG. 3: (a) Magnetic field dependence of the spin polarization parameter S at the equilibrium state 

for the cases T = MeV (solid curve), T =10 MeV (dashed curve) and T = 20 MeV (dashdot 
curve), and a fixed value of the density, p = 0.2 fm~^. (b) Same as in the top panel but for the 
ratio T/e*f. 
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FIG. 4: Free energy per particle at the equilibrium state as a function of: (a) the density p for 

the cases T = MeV (solid curve), T = 10 MeV (dashed curve) and T = 20 MeV (dashdot 

curve), and a fixed value of the magnetic field, B = 5x 10^* G, (b) the temperature T for the cases 

B = G (solid curve), B = 10^^ G (dashed curve) and B = 5 x 10^^ G (dashdot curve), and a 

fixed value of the density, p = 0.2 fm"^, (c) the magnetic field B for the cases T = MeV (solid 

curve), T = 10 MeV (dashed curve) and T = 20 MeV (dashdot curve), and a fixed value of the 
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density, p = 0.2 fm~^. 
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FIG. 5: Pressure of spin polarized neutron matter as a function of: (a) the density p for the cases 

B = 0G (sohd curve), B = 10^^ G (dashed curve) and B = 5x 10^^ G (dashdot curve), and a fixed 

value of the temperature, T = 20 AleV, (b) the density p for the cases T = MeV (solid curve), 

r = 10 MeV (dashed curve) and T = 20 MeV (dashdot curve) , and a fixed value of the magnetic 

field, B = 5 X 10^^ G, (c) the temperature T for the cases B = G (solid curve), B = 10^^ G 

(dashed curve) and B = 5x 10^^ G (dashdot curve), and a fixed value of the density, p = 0.2 fm"^. 
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FIG. 6: Magnetic field dependence of tlie effective masses of spin-up (solid curve) and spin-down 
(dashed curve) neutrons corresponding with the equilibrium state at the fixed values of the tem- 
perature, T = 10 MeV, and the density, p = 0.3 fm~^. 
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